Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

A physical decomposition of the stress tensor for complex flows

Traditionally, the components of the stress with respect to a relevant coordinate system are used for the purpose of stress visualisation and interpretation. A case for using a flow dependent measure to interpret and visualise stress is made for two dimensional flow, together with a suggestion for extending the idea to three dimensions. The method is illustrated for Newtonian and Oldroyd B fluids in both the eccentrically rotating cylinder and flow past a cylinder benchmark problems. In the context of a generalised Newtonian fluid, the relation between the flow-dependent stress measure to other field variables under certain flow conditions, is examined and is indicative of its importance in complex flow.

Integral operators commuting with dilations and rotations in generalized Morrey-type spaces

We find conditions for the boundedness of integral operators $K$ commuting with dilations and rotations in a local generalized Morrey space. We also show that under the same conditions, these operators preserve the subspace of such Morrey space, known as vanishing Morrey space. We also give necessary conditions for the boundedness when the kernel is non-negative. In the case of classical Morrey spaces, the obtained sufficient and necessary conditions coincide with each other. In the one-dimensional case, we also obtain similar results for global Morrey spaces. In the case of radial kernels, we also obtain stronger estimates of $Kf$ via spherical means of $f$. We demonstrate the efficiency of the obtained conditions for a variety of examples such as weighted Hardy operators, weighted Hilbert operator, their multidimensional versions, and others.     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象，发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式，避免在单元内进行速度平均，从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比，网格的边界通量由节点处的多维黎曼求解器提供，克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增，给出稀疏波出现的合理判据，从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象，同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.     

A relaxation Riemann solver for compressible two-phase flow with phase transition and surface tension

The dynamics of two-phase flows depend crucially on interfacial effects like surface tension and phase transition. A numerical method for compressible inviscid flows is proposed that accounts in particular for these two effects. The approach relies on the solution of Riemann-like problems across the interface that separates the liquid and the vapour phase. Since the analytical solutions of the Riemann problems are only known in particular cases an approximative Riemann solver for arbitrary settings is constructed. The approximative solutions rely on the relaxation technique.The local well-posedness of the approximative solver is proven. Finally we present numerical experiments for radially symmetric configurations that underline the reliability and efficiency of the numerical scheme.     

Cartan Invariants in the Category of Restricted Representation for gl(m|n)

Let k be an algebraically closed field of characteristic p 2, and gl(m|n) be the general linear Lie superalgebra over k. The Cartan invariants for the restricted supermodule category for gl(m|n) are presented.     

An embedded boundary framework for compressible turbulent flow and fluid–structure computations on structured and unstructured grids

The finite volume method with exact two‐phase Riemann problems (FIVER) is a two‐faceted computational method for compressible multi‐material (fluid–fluid, fluid–structure, and multi‐fluid–structure) problems characterized by large density jumps, and/or highly nonlinear structural motions and deformations. For compressible multi‐phase flow problems, FIVER is a Godunov‐type discretization scheme characterized by the construction and solution at the material interfaces of local, exact, two‐phase Riemann problems. For compressible fluid–structure interaction (FSI) problems, it is an embedded boundary method for computational fluid dynamics (CFD) capable of handling large structural deformations and topological changes. Originally developed for inviscid multi‐material computations on nonbody‐fitted structured and unstructured grids, FIVER is extended in this paper to laminar and turbulent viscous flow and FSI problems. To this effect, it is equipped with carefully designed extrapolation schemes for populating the ghost fluid values needed for the construction, in the vicinity of the fluid–structure interface, of second‐order spatial approximations of the viscous fluxes and source terms associated with Reynolds averaged Navier–Stokes (RANS)‐based turbulence models and large eddy simulation (LES). Two support algorithms, which pertain to the application of any embedded boundary method for CFD to the robust, accurate, and fast solution of FSI problems, are also presented in this paper. The first one focuses on the fast computation of the time‐dependent distance to the wall because it is required by many RANS‐based turbulence models. The second algorithm addresses the robust and accurate computation of the flow‐induced forces and moments on embedded discrete surfaces, and their finite element representations when these surfaces are flexible. Equipped with these two auxiliary algorithms, the extension of FIVER to viscous flow and FSI problems is first verified with the LES of a turbulent flow past an immobile prolate spheroid, and the computation of a series of unsteady laminar flows past two counter‐rotating cylinders. Then, its potential for the solution of complex, turbulent, and flexible FSI problems is also demonstrated with the simulation, using the Spalart–Allmaras turbulence model, of the vertical tail buffeting of an F/A‐18 aircraft configuration and the comparison of the obtained numerical results with flight test data. Copyright © 2014 John Wiley & Sons, Ltd.     

Strategies to cure numerical shock instability in the HLLEM Riemann solver

The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver that is known to be plagued by various forms of numerical shock instability. In this paper, we clarify that the shock instability exhibited by this scheme is primarily triggered by the spurious activation of the antidiffusive terms present in the first and third Riemann flux components on the transverse interfaces adjoining the shock front due to numerical perturbations. These erroneously activated terms are shown to counteract the favorable damping mechanism provided by its inherent HLL-type diffusive terms, causing an unphysical variation of the conserved quantity ρu both along and across the numerical shock. To prevent this, two distinct strategies are proposed termed as S elective W ave M odification and A nti D iffusion C ontrol. The former focuses on enhancing the quantity of the favorable HLL-type dissipation available on these critical flux components by carefully increasing the magnitudes of certain nonlinear wave speed estimates, while the latter focuses on directly controlling the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and to estimate a von Neumann–type stability bounds on the CFL number associated with their use. Results from a variety of classic shock instability test cases show that the proposed strategies are able to provide excellent shock stable solutions even on grids that are highly elongated across the shock front without compromising the accuracy on inviscid contact or shear dominated viscous flows.     

Conservation Laws of The Generalized Riemann Equations

Two special classes of conserved densities involving arbitrary smooth functions are explicitly formulated for the generalized Riemann equation at arbitrary N. The particular case when N = 2 covers most of the known conserved densities of the equation, and the result is also extended to the famous Gurevich-Zybin, Monge-Ampere and two-component Hunter-Saxton equations by considering certain reductions.